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Abstract. In recent years, three exceptional discretizations of the </> 4 theory have 
been discovered |SW94I IBT97I lK03 which support translationally invariant kinks, i.e. 
families of stationary kinks centred at arbitrary points between the lattice sites. It 
has been suggested that the translationally invariant stationary kinks may persist as 
sliding kinks, i.e. discrete kinks travelling at nonzero velocities without experiencing 
any radiation damping. The purpose of this study is to check whether this is indeed the 
case. By computing the Stokes constants in beyond-all-order asymptotic expansions, 
we prove that the three exceptional discretizations do not support sliding kinks for 
most values of the velocity — just like the standard, one-site, discretization. There are, 
however, isolated values of velocity for which radiationless kink propagation becomes 
possible. There is one such value for the discretization of |SW94| and three sliding 
velocities for the model of IKf 



PACS numbers: 05.45.Yv, 63.20.Pw 
1. Introduction 

Spatially discretized partial differential equations (or, equivalently, chains of coupled 
ordinary differential equations) have attracted considerable attention recently. One of 
the issues that has been vigorously debated and that will concern us in this paper, is 
whether discrete systems can support solitary waves travelling without losing energy to 
resonant radiation and decelerating as a result. We address this issue for one of the 
prototype models of nonlinear physics, the 4 -theory: 

u tt = u xx + ^u(l-u 2 ). (1.1) 
The </> 4 -equation (jl.lj) is Lorentz- invariant, and so the existence of the travelling kink 

u(x, £) = tanh- — t , (1.2) 
*{ On sabbatical leave from the University of Cape Town, South Africa 
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where |c| < 1 and s 6 1, is an immediate consequence of the existence of the stationary 
kink for c = 0. On the other hand, if we discretize equation (jl.lj) in x, 

u n+l — 2tt n + v , v 

= ?2 h/(Wn-l,W n ,W n+ l), (1.3) 

the translation and Lorentz invariances are lost, and the existence of the travelling kink 
(and even of an arbitrarily centred stationary one) becomes a nontrivial matter. In 
equation 1)1.3)1 . u n G M, n G Z, t G 1, /i is the lattice spacing, and the nonlinearity 
J(it n _i, u n , u n+ i) satisfies the continuity condition 

/(«,«,«) = \u{l-u 2 ). (1.4) 
We restrict ourselves to symmetric discretizations, i.e. 

/(«n-l, «n, M n+ l) = f( u n+l, U m U n -i). (1.5) 

Equation ()1.1)1 results from ()1.3)1 in the continuum limit, where u n (t) = u(x n , t), x n = nh 
and h — > 0. In this limit, the truncation error of the Taylor series is 0(h 2 ). We shall 
be concerned with monotonic kink solutions of ()1.3)) : u n+ i(t) > u n (t) for all n € Z. As 
/i — >• 0, such monotonic discrete kinks approach the continuous kink ([1.2)1 . 

The most common, one-site discretization of the nonlinearity function is given by 

f(u n -!, u n , u n+ i) = -u n {\ - u 2 n ). (1.6) 

It is a well established fact )L88j . however, that the discrete Klein-Gordon equation 
([1.3)1 + ([1.6)) admits only a countable set of stationary monotonic kinks with the boundary 
conditions 

lim u n {t) = —1, lim u n {t) = +1. (1.7) 

n— >— oo n— >+oo 

Physically, this fact is related to the presence of the Peierls-Nabarro barrier, an effective 
potential periodic with the spacing of the lattice. Half of the stationary kinks are centred 
at the minima (the on-site kinks) and the other half (the off-site kinks) at the maxima 
of the Peierls-Nabarro potential. There are no continuous families of stationary discrete 
kinks of the form u n = u(n — s), with s a free parameter, which would interpolate 
between the two solutions. In an abuse of terminology, we will be calling such families 
"translationally invariant kinks" — although, in the first place, translation invariance 
is a property of an equation rather than a solution, and in the second, all lattice 
equations are of course not translationally invariant. As for propagating waves, of special 
importance are kinks moving at constant speed and without the emission of radiation. 
We will be referring to such kinks, i.e. solutions of the form u n = u(n — ct — s) where 
w(£) is a monotonically growing function satisfying the boundary conditions (Jl. 7)1 . as 
sliding kinks, to emphasise the fact that they do not experience any radiative friction. 
Being an obstacle to the "translationally invariant" kinks, the Peierls-Nabarro barrier 
is also detrimental to the existence of sliding kinks — at least for small c (see reviews 
in [SB] and $W\ ). 
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In an attempt to find a discrete model with "translationally invariant" and sliding 
kinks, Speight and Ward |SW941 IS97j considered a hamiltonian discretization of the 



form 



f{u n - 1 ,u n ,u n+1 ) = -(2« t! + M n+1 ) II 



+ — (2w n + w n _i) I 1 I. (1.8) 

In the static limit, the corresponding energy admits a topological lower bound which is 
saturated by a first- (rather than second-) order difference equation. This equation is 
readily shown to have a one-parameter continuous family of stationary kink solutions 
u n = u(n — s) for < h < 2 (see Proposition 1 in |S97p . The parameter s of the family 
defines the position of the kink relative to the lattice. Since all members of the family 
have the same (lowest attainable) energy, the stationary kink experiences no Peierls- 
Nabarro barrier. As for travelling kinks, Speight and Ward's numerical simulations 
revealed that although moving kinks in this model do lose energy to Cherenkov radiation 
and decelerate as a result, this happens at a slower rate than a similar process in equation 
(ll.fij) (see figures 4 and 5 in |S97j ). 

Another line of attack was chosen by Bender and Tovbis |BT97j who proposed a 
different discretization supporting a continuous family of arbitrarily centred stationary 
kinks: 

/(n n _i, u n , u n+1 ) = - (u n+ i + u n _i) (l-u 2 n ) . (1.9) 

In this case, the family arises due to the suppression of the stationary kink's resonant 
radiation. In fact, the family of stationary kinks can be found explicitly as 

u n (t) = tanh [a{n — s)] , (1-10) 

where a = arcsinh (h/2) for all /iGl. (The solution (|1.10|) coincides with the stationary 
dark soliton of the repulsive Ablowitz-Ladik equation |HA93j .) 
Finally, the nonlinearity 

f(u n -!, U n , U n+ i) = - (u n+1 + U n -i) (2 - u 2 n+1 - • (1-11) 



was introduced by Kevrekidis |K03j . who demonstrated the existence of a two-point 
invariant and hence a first-order difference equation associated with the stationary 
equation. Consequently, the discretization (jl.llj) also supports a continuous family 
of stationary kinks for all h G [0, ho] with some ho > 0. (For general discussion, see 
[S99], |B()P05j and |DKY05a] .) A relevant property of the model fll.llj) . which is related 
to the existence of a two-point invariant |K03j and indicates some additional underlying 
symmetry, is the conservation of momentum. (See also |DKY05b] .) 

Since the reasons for the nonexistence of "translationally invariant" kinks and 
of sliding kinks are apparently related (the breaking of symmetries of the underlying 
continuum theory or, speaking physically, the presence of the Peierls-Nabarro barrier), 
the availability of "translation- invariant" stationary kinks in the models ()1.8|) . ()1.9|) and 
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(II. lip suggests that they might have sliding kinks as well. It is the purpose of the present 
study to find out whether this is indeed the case. We shall analyse the persistence of 
continuous families of stationary kinks u n = u(n — s) for nonzero velocities; in other 
words, examine the existence of solutions of the form u{n — ct — s) where u(z) is a 
monotonically growing function satisfying (|1.7jl . and c ^ 0. We develop an accurate 
numerical test in the limit h — > which shows whether or not standing and travelling 
kinks of the discrete 4 model (jl.3|) bifurcate from the exact kink solutions (jl.2|) of its 
continuous counterpart (jl.lj) . The analysis of this bifurcation poses a singular problem 
in perturbation theory which can be analysed using two (inner and outer) matched 
asymptotic scales on the complex plane |TT.I98) iTOOaj . In particular, the nonvanishing 
of the Stokes constant in the inner asymptotic equation serves as a sufficient condition 
for the non-existence of continuous solutions of the difference equations |TT.T98j . 

Our test will be based on computing the Stokes constant for the differential- 
difference equation underlying the lattice system. We will examine all four 
discretizations of the 4 theory mentioned above, i.e. equations (jl.fi)) . (jl.8j) . (jl.9j) and 
(11.11)1 . Since translationally invariant stationary kinks u n = u(n — s) do exist for the 
three exceptional nonlinearities ()1.8j) . (jl.9)l and (jl.llj) . the Stokes constant is a priori 
vanishing for c = in these three cases. However, we will show that in all three 
cases the Stokes constant acquires a nonzero value as soon as c deviates from zero. 
It remains nonzero for all c except a few isolated values which define the particular 
velocities of the sliding kinks in the corresponding model. There is one such isolated 
zero of the Stokes constant for the nonlinearity (jl.8j) and three sliding velocites for the 
discretization (|1.11)) . Consequently, the main conclusion of this work is that the sliding 
kinks, i.e. kinks travelling at a constant speed without the emission of radiation, can 
occur only at particular values of the velocity. The sliding velocities are, of course, 
functions of the discretization spacing h, so that sliding kinks arise along continuous 
curves on the (c, n)-plane. 

We conclude this introduction with a remark on a convention adopted in the 
remainder of this paper — namely, that the linear part of the function f(u n -i, u n , u n+ i) 
in (J1.3)) can always be fixed to \u n without loss of generality. Indeed, the most general 
function satisfying (J1.4j) and (jl.5j) is / = (| — 2a) u n + a (u n+1 + u n _i) + cubic terms, 
where a is arbitrary. Since h 2 in (jl.3j) is also a free parameter, we can always make a 
replacement h — > h such that 1/h 2 + a = 1/h 2 . This gives 

/(n n _i, u n , u n+ i) = -u n - Q{u n -!, u n , u n+1 ), (1.12) 

where Q is a homogeneous polynomial of degree 3 which is independent of the parameter 
h. 

The plan of this paper is as follows. In the next section (section EJ) we review the 
construction of the outer and inner asymptotic solutions in the limit h — > 0. Section 
El contains details of the numerical computation of the Stokes constants while the last 
section (section 0]) summarises the results of our work. 
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2. Inner and outer asymptotic expansions in the limit h — > 

We are looking for a sliding-kink solution of the discrete 4 models (jl.Hj) in the form 

u n (t) = <f)(z), z — h(n — s) — ct, (2.1) 

where <fi(z) is assumed to be a twice differentiable function of z £ JR., that satisfies the 
differential advance-delay equation 

c (2) = — + -0(z) - Q {(j){z - h), (p(z), (j){z + h)) , (2.2) 

with the boundary conditions <fi(z) — > ±1 as z — > ±00. The velocity c is assumed 
to be smaller than 1 in modulus. If a solution to this boundary-value problem (i.e. 
a heteroclinic orbit) exists, then the parameter s is arbitrary due to the translation 
invariance of the advance-delay equation ()2.2j) . The scaling parameter h (which stands 
for the lattice step-size) can be used to reduce equation (J2.2j) to a singularly perturbed 
differential equation as h — > |TTJ98j . Formal asymptotic solutions of the problem 
(12. 2 j) can be constructed at the inner and outer asymptotic scales. The formal series 
represent convergent asymptotic solutions of the singular perturbation problem only if 
the Stokes constants are all zero |T00a| . 

Asymptotic analysis beyond all orders of perturbation theory was pioneered by 
Kruskal and Segur |KS91j and has been utilised by many authors. It was extended by 
Pomeau et. al. |PRG88| to allow the computation of radiation coefficients from Borel 
summation of series rather than from the numerical solution of differential equations. 
Essentially the same method has been applied to different problems by Grimshaw and 
Joshi |HJ95l IH95] and Tovbis and collaborators |TT,)981 rTUUal iTnnbl ITP05| . In this 
paper, we shall work with formal inner and outer asymptotic series for the problem ()2.2j) 
without attempting rigorous analysis of their asymptoticity. 

2.1. Outer asymptotic series 

Assuming that the solution <p(z) is a real analytic function of z, we consider the Taylor 
series for the second difference in a strip V$ = {z £ C : |Imz| < 5}, where 5 > 0: 

00 g 

<j>{z + h)- 2#z) + <j>{z -h) = h 2 <f>"(z) + g h 2n -^^ n \z). (2.3) 

Since the cubic polynomial Q(u n -i, u n , u n+ i) satisfies the continuity and symmetry 
relations (jl.4j) and ()1.5|) . the nonlinearity of (J2.2j) can also be expanded in a Taylor 
series in the same strip: 

^ 00 

Q (<P(z - h), 0(z), fa + h)) = -<p 3 (z) + h 2n Q2n (0, (<P') 2 , {2n) ) , (2.4) 

n=l 

where the coefficients Q2 n depend on even derivatives and even powers of odd derivatives 
of 4>(z), and also Q2n(0, 0, ...,0) = 0. The differential advance-delay equation ()2.2|) can 



Travelling kinks in discrete <p A models 



6 



thus be written as 

(1 - <?W + ^(1 -0 2 ) + E ( (2^2), ^ (2n+2) - (<^> (^') 2 , ■ ■ • , (2n) )) = 0. 

(2.5) 

For h = 0, equation (|2.5|) becomes the travelling wave reduction of the continuous model 
(II. lj) . with the explicit solution 

(f> (z) = tanh£; e = - 7^== , |c| < 1. (2.6) 

2vl — 

We will search for solutions of equation (|2.5jl of the form 

oo 
n=l 

Substituting the expansion (|2.7|) into (|2.5|) we get, at order /i 2n , 

£>4>2n = H2ni 

where the linearised operator £ is given by 

£ = -^ + 4-6sech 2 £, 

and H 2n are polynomials in <fr , <p 2 , ■ ■ ■ , <p2n-2 and their derivatives. The kernel of C is 
one-dimensional, and spanned by an even eigenfunction yo = sech 2 £. The rest of the 
spectrum of C is positive. It is not difficult to prove by induction that if 4>2k{z) are 
all odd in z for k = 0, 1, ...,n — 1, the nonhomogeneous term H 2n is also odd in z and 
hence, by the Fredholm alternative, there exists a unique odd bounded solution 4>2n{z) 
for z G R. Moreover, since i?2n decays to zero exponentially fast as \z\ — > oo, the 
function <p2n{z) is also exponentially decaying for any n > 1. The perturbation 02 (2) in 
particular satisfies the nonhomogeneous equation 

£02 = (</£ v) + o^o' + + 70o(0o) 2 ) , (2.8) 

where the numerical coefficients depend on whether the nonlinearity function / is given 
by (USD, (HID, (USD or (HID): 

One-site (JUBJ) : a = /3 = 7 = 0; 

Speight- Ward (EED : a = 1, /3 = 7 = -4; 

Bender-Tovbis ((131) : a = 3, /3 = -3, 7 = 0; 

Kevrekidis (fl~Xl"| : a = 3, /? = -9, 7 = -6. 

The odd bounded solution 2 {z) of the nonhomogeneous equation ()2.8|) is: 

2 (z) = .4tanh£sech 2 £ + 23£sech 2 £, (2.9) 

where 

(l-c 2 )( 7 + 2/j) + 6 (l-c 2 )(a + /3) + l 

72(1 -c 2 ) 2 ' 24(1 -c 2 ) 2 

The hat in the series (|2.7|) indicates that the series is formal, i.e. it may or may not 
converge |TTJ98| fTOOaj . depending on the choice of c and Q in equation (J2.2|) . We shall 



be referring to (|2.7|) as the outer asymptotic expansion. 



Travelling kinks in discrete 4 models 7 
2.2. Inner asymptotic series 

The leading-order term ()2.6|) of the outer expansion (|2.5J) has poles at £ = y(l + 2n), 
where fiGZ. We apply the scaling transformation 

z = h( + inVl - c 2 , 0(z) = -^(C) (2.10) 

to equation ()2.2|) in order to study the convergence of the formal asymptotic solution 
(1271) near the pole f = f (see jTT,)98L ITOOaj ). This yields the following differential 
advance-delay equation for 

c y (0 = ^(c + 1) - 2^(o + v(c - 1) - q m - 1), vco, v>(c + 1)) + |W (2.11) 

The following are the cubic functions Q for each of the four discretizations that we deal 
with in this paper: 

One-site (H2J) : Q = ^ 3 {(); 

Speight- Ward (HD : Q = [^ 3 (C + 1) + 3^ 2 (C + l)V(C) + 3^(C + 1)^ 2 (C) 

+4^ 3 (0 + 3^(C - i)^ 2 (C) + 3^ 2 (C - l)^(C) + ^ 3 (C - i)] ; 

Bender-Tovbis flH| : Q = ^ 2 (C) bP(C + 1) + <KC ~ 1)] I 

Kevrekidis (PU : Q = - [</> 3 (C + 1) + V> 2 (C + 1)V(C ~ 1) + V>(C + 1)^ 2 (C - 1) 

+^ 3 (C - 1)] • 

We note that the heteroclinic orbit becomes small as h — > under the normalization 
(|2.1()jl : if <f>(z) — > ±1 as z — > ±oo, then V(C) ~~ ^ ^ as Re£ ^ ±oo. The formal 
asymptotic series ()2.7|) in the new variables (|2.10j) becomes a new formal series 

oo 
n=l 

where each term V\i(C) can be expanded in a formal series in descending powers of £■ 
In particular, the leading-order function "0o(O nas the general form 

oo 

«0 = E^;. (2-13) 

n=0 S 

By comparing the series (J2.12|) and (J2.13j) with the solutions (J2.6p and (J2.9|) in variables 
(I2.10|) . we note the correspondence: 

„ , pr—s (l-e 2 )( 7 + 2,9)+6 

We shall be referring to ()2.12|) as the inner asymptotic expansion. The odd powers of 
h in the inner asymptotic expansion ()2.12j) appear due to the matching conditions with 
the outer asymptotic expansion ()2.7|) under the scaling (|2.10jh as well as due to the 
non-zero boundary conditions for the heteroclinic orbits ip{Q - > ±^ as Re£ — > ±oo. 
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2.3. Leading-order problem for an inner solution 

Convergence of the formal inverse-power series ()2.13)1 for the leading-order solution 
^o(C) depends on the values of the Stokes constants |T00a| . Computation of the Stokes 
constants is based on Borel-Laplace transforms of the inner equation (gTTJ) |TT,)98j . 
Assuming continuity in h, we study the leading-order solution ^o(C) = h m /i->o V'(C) °f 
the truncated inner equation 

c 2 <(C) = MC + 1) - 2MC) + MC - 1) - Q (MC - 1), MC),MC + 1)) • (2.14) 

By substituting the series (j2.13j) into equation (|2.14jl . one can derive a recurrence relation 
between the coefficients in the set {a n }^L . The Stokes constants can be computed from 
the asymptotic behavior of the coefficients a n for large n. Alternatively, the leading- 
order solution "00 (C) an d the Stokes constants can be defined by using the Borel-Laplace 
transform: 

MC) = / V (p)e-^dp. (2.15) 



The choice of the contour of integration 7 determines the domain of ipo(C) i n t ne complex 
C-plane. We define two solutions and V'o (Oj which lie on the stable and unstable 

manifolds respectively, such that 

lim #(()=0, lim 4 U) (O = 0. (2.16) 

Kei,— >+oo Ref— >— 00 

We note that the stable and unstable solutions tend to the stationary point at the origin, 
since the heteroclinic orbits connect the stationary points at ip — ±h which move to the 
origin as h — > 0. The three stationary points coalesce to become a degenerate stationary 
point at the origin within the truncated inner equation 1)2.14)1 . 

The Borel-Laplace transform ()2.15)) produces the stable solution i/Jo\0 when the 
contour of integration 7^ lies in the first quadrant of the complex p-plane and extends 
from p = to p = 00. Similarly, it produces the unstable solution ip^iC) when the 
contour of integration 7 n lies in the second quadrant. We choose the integration contours 
in such a way that argp — > tt/2 as p — > 00, so that the solutions ^0 (0 an d ^0 (C) are 
defined by ()2.15)) for all complex ( with Im£ < 0. 

The Borel transform Vo(p) satisfies the following integral equation, which follows 
from flZ31D and 

(4 sinh 2 \ - cV) Vb(p) = Q [Vo(p)] . (2.17) 

Here, Q [V(p)} denotes a double convolution of V(p) with itself (in this case, the hat is 
used to denote an operator). We list below the convolutions Q [V(p)] for each of the 
four models under consideration: 

One-site : 2Q = V(p) * V(p) *V(p); 

Speight-Ward dUD : 18Q = coshp [V(p) * V(p) * V{p)\ + 3 {coshp [V(p) * V{p)}} 

* V{p) + 3 [coshp V(p)\ * V(p) * V(p) + 2V(p) * V(p) * V{p); 
Bender-Tovbis dUJ) : 2Q = [coshp V(p)} * V{p) * V(p); 
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Kevrekidis (ITTTT1) : 4Q = coshp [V(p) * V(p) * V{p)} 

+ [coshp V(p)\ * [e p V(p)} * [e~ p V(p)] , 
where the asterisk * denotes the convolution integral for the Borel-Laplace transform: 

v(p) * w(p) = r vip - P1 )w( P1 )d P1 , 

Jo 

and the integration is performed from the origin to the point p on the complex plane, 
along the contour 7. The inverse power series ()2.13|) for the limiting solution ^o(C) 

becomes the following power series for the Borel transform Vo(p): 

00 

Vob) = 5>n/ n , ^ = ?flp ( 2 - 18 ) 

n=0 ^ ' ' 

where Vq = a = 2yi — c 2 . The hat denotes a formal series which might only converge 
for some values of p. The virtue of the integral form (|2.17|) is that the limiting behavior 
of v n for large n can be related to singularities of Vq( P ), which in turn correspond to the 
oscillatory tails of "00(C)- 

If the sliding kink exists, the inverse-power series for ipo(C) w iU converge for all 
( 6 C such that Im£ < 0. This implies that the stable and unstable solutions "0^(0 
and ^o^(C) coincide, i.e. that the contour 7^ in the right half of the complex p-plane 
can be continously deformed to the contour 7^ in the left half-plane (see figure QJ. 
If, however, there are any singularities between the two contours, then a continuous 
deformation is possible only if the residues are zero. The residues are proportional to 
the values of the Stokes constants. When the Stokes constants are nonzero, the formal 
power series (j2.18J) for the solution V (p) of the integral equation (|2.17|) diverges for 
some values of p in the sector between the contours 7^ and j u . 
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2.4. Stokes constants 

The Borel transform Vo(p) is singular near the points in the p-plane where the coefficient 
in front of Vo(p) on the left-hand side of the integral equation (|2.17j) vanishes |T00aj . 
except for the point p = where the right hand side is also zero. That is, singularities 
occur when (2/p) sinh(p/2) = ±c. The location of these singularities is important 
because the stable and unstable solutions are not, in fact, uniquely defined by ([2.16)1 : 
different solutions are generated depending on where the contours lie relative to the 
singularities of Vq(p) with Rep 7^ 0. Exploiting this nonuniqueness, we wish to choose 
the contours 7 S and 7^ to lie above all the singularities with nonzero real part; this will 
minimise the number of singularities between the stable and unstable solutions. 

It is not difficult to show that the contour 7,, extending from to 00 can be chosen 
in such a way, i.e. so that that there are no singularities between it and the imaginary 
axis. Indeed, assume, for definiteness, that c > 0. Let (2n c — 1) be the number of 
positive roots of the equation sing = cq and denote the real and imaginary parts of p/2 
by k and q: p/2 = k + iq. In the (k, g)-plane, consider a rectangular region T> bounded 
by the horizontal segments q = 2%n and q = e at the top and bottom, and vertical 
segments k = — e and k = e on the left and right. Here n is any positive integer greater 
than n c and e > is taken to be small. Using the argument principle, we can count the 
number of (complex) roots of the equation sinh(p/2) = c(p/2) in the region T>. We have 

p cosh k sin q — cq 

tan arg - = — . 

2 sinn k cos q — ck 

On the right lateral side, where k — e, this becomes 

p 1 sing -cq 

tan arg - « . (2.19) 

2 e cos q — c 

As we move from q = e to q = 27m, the numerator in ()2.19|) will change sign (2n c — 1) 
times. In a similar way, moving down along the left side there will be (2n c — 1) more 
zero crossings, while no zero crossings will occur along the horizontal segments. This 
means that the argument can change by no more than (4n c — 2)ir and hence there are 
at most (2n c — 1) roots in the region D, no matter how large n is. Similarly, we can 
show that the equation sinh(jo/2) = — c(p/2) has no more than 2n c roots in the region 
T>, if 2n c is the number of positive roots of sing = —eg. The upshot is that for any finite 
c, there are only a finite number of singularities with small real parts; the singularities 
cannot accumulate to the imaginary axis. For c ^ 0, the singularities with nonzero real 
parts lie on the curves 



1 . , „ ,0 . 1 



q = ±\ — cosh 2 k — k 2 coth 2 n — ► ±- cosh k as \k\ 



c 2 c 



00. 



Accordingly, in order for the integration contours 7 S and 7„ to lie above these 
singularities, they must be curvilinear (and not just rays) as shown in figure ^ 

Having chosen the contours 7^ and 7„ to lie above the singularities in the first and 
second quadrants respectively, the only singularities of V (p) that determine whether 
the stable solution V'o^(C) can be continuously transformed into the unstable solution 
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ijj^ (() are those at non-zero pure imaginary values of p. We will be referring to these 
values as resonances. The set of resonances TZ C is defined by the transcendental equation 

TZ C = jp = ik, k G R+ : | sin ^ = ±c| . (2.20) 

When c = 0, the set TZo is infinite-dimensional and can be described explicitly: 

TZ = {p = 2vmi, n£N}. 

Let p\ = iki be the smallest imaginary root in the set TZ C . It is clear from ()2.20j) that 
< k\ < 2tt for c G (0, 1), so that k\ — > 27r as c — > + and fci — > as c — > 1~. The 
set of resonances TZ C is finite-dimensional for c G (0, 1) and it consists of only one root 
Pi = iki for c G (ci, 1), where ci ~ 0.22. 

Due to the resonances, a function ^o(C) that satisfies the truncated inner equation 
(12.14)1 may have oscillatory tails as [Re Cl ~~ > °°- Adding the solutions of equation (J2.14jl 
linearised about ^o(C); the general bounded solution of (j2.14j) in the limit |C| — ^ oo can 
be represented as |TTJ98j : 

^ (() = ^ (£) + ^ a n n (()e~ ikn<: + multiple harmonics. (2.21) 

ik n eTZ c 

Here, V>o(C) is given by the power series (|2.13jl : a n are coefficients which we will be 
referring to as amplitudes in what follows; k n > are roots of (J2.20jl for p = ik n , and 
the functions (p n (()e~ lkn< * , n > 1, satisfy the linearised truncated inner equation ([2.14)1 . 
In particular, the equation for the leading-order term <fi(C) is 

e" ifc Vi(C + 1) + {c 2 k\ - 2)^(0 + e ife Vi(C - 1) + 2ic 2 k 1 0' 1 (C) - cVi(C) 

= DiQ yMC - 1) + D 2 Q ^(C) + D 3 Q X (C + 1), (2.22) 

where Di 2 ^Q are the partial derivatives of Q(^o(C ~~ 1)> V'o(O) ^o(C + 1)) with respect 
to its first, second and third argument respectively, evaluated at ^ = ^o(C)- 

If the amplitude a n is nonzero for some n, the formal power series ()2.13|) does not 
converge because the solution (J2.21)l does not decay as |Re£| — > oo. The amplitudes a n 
are proportional to the Stokes constants computed for the formal power series (J2.13)) . 
Each oscillatory term in the sum ([2.21J1 becomes exponentially small in h when we 
transform from ( to z using the transformation (|2.1()j) . Since p\ = iki is the element of 
1Z C with the smallest imaginary part, it follows that the n — 1 term dominates the sum 
in (|2.21j) when the transformation ()2.10|) is made (unless a\ — 0). Furthermore, when 
c G (ci, 1), where C\ ~ 0.22, it is the only term in the sum since the resonant set 1Z C 
consists of just the one root p 1 = iki. We shall, therefore, only be concerned with the 
leading-order Stokes constant, which multiplies the function <^i(C). 

If "00 (C) is given by the power series (|2.13|) . the solution of the linearized equation 
(12.22)1 can also be represented by a formal power series: 



oo 



M0 = C}M-\ (2.23) 
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where we can set b = 1 due to the linearity of (|2.22|) . Substituting ([2.13)1 and ([2.2 3)) 
into ()2.22|) and using (J2.20)) . the coefficient bi can be determined from 



2irC r_1 (c 2 h - sinh) + C~ 2 [r(r - l)(cosA;i - c 2 ) 

+2i6i(r - l)(c 2 fci - sin fci) - 6(1 - c 2 )] + C(C 



0. 



In this equation, the coefficient of each power of ( should be set to zero. In order to set 
the coefhcent in front of the first term to zero in the situation where c ^ 0, we must 
choose r = 0. The second term then gives 



bi 



3i(l 



c 2 k\ — sin k\ 



after which all the other coefficients 62, ^3, • • • , can be computed recursively. On the 
other hand, in the situation with c = 0, we have ki = 2n and the coefficient in front 
of C, r ~ l is zero regardless of the value of r. Setting the coefficient in front of ( r ~ 2 to 
zero requires that we choose either r = 3orr = — 2, and hence we have two different 
descending-power series, one starting with £ 3 and the other one with (~ 2 . We shall focus 
on the former as it dominates the latter in the limit ( — ► 00. Again, the succeeding 
terms in ()2.23|) are determined recursively. 

Thus, we have established that the leading-order oscillatory term in the expansion 
(12.21)1 behaves as 



-ifciC 



for c 7^ and 
for c = 0. 



(2.24) 



For c 7^ 0, the two leading order terms in the expression above are generated by, 
respectively, a simple pole and a logarithmic singularity of the Borel transform Vo(p) at 
p — iki. For c = they are generated by a quadruple pole of Vo(p) at p = 2ni. From 
the fact that Vq is an even function of p, we deduce the structure of this function near 
the poles: 

' klK^c) 



Vo(p) 



p 2 + k 2 -°(c)Hp 2 + kl) + 



6(27r) 8 ^ 



2(2vr) 6 p 



+ . . . 



(for c^0) 
(for c = 0) 



(2.25) 



(p2 + 47r 2)4 ^ 2 + 4?r 2)3 

as p — > ±ifci. Here Ki(c) and Si are the leading-order Stokes constants for c ^ and 
c = 0, respectively; cr(c) and p are independent of p, and . . . stands for terms with even 
slower growth as p — >• ik\. 

To show that these singularities do indeed give rise to the oscillatory tails in 1)2.24)1 . 
we compare the two integrals ip^ and for a given value of (. To this end, we 
deform the paths of integration 7 S and 7 n to 7^ and 7^ respectively, without crossing 
any singularities. This is illustrated in figure ^ There are two contributions to the 
difference V'o (C) — V'o (0- The first comes from integrating around the pole at p — iki, 



y o is; s^o 

and is equal to 27ri times the residue of the function Vo(p)e~ p< * at p 



ik\, determined 
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from (|2.25|) . The second contribution (manifest only in the c 7^ case) arises because 
the integrand increases as the singularity is encircled, since it is a branch point of the 
logarithm. Since In z can be written as In \z\ + iargz, where z = p — iki, we see that 
Vq(p) increases by — 2iria(c) as the branch point p = ik\ is encircled in the c / 
case. Therefore, the difference in the integrand of (|2.15jl along the paths 7^ and 7^ is 
— 2iria(c)e~ p< * , which must be integrated along the path of integration from p = ik\ to 
infinity, to give — 27ricr(c)e _lfcl ^/C. (We have considered the integration on a Riemann 
surface in order to account for branch points.) Adding together the two contributions 
discussed above, we have 

2nia(c) 



7rkiKi(c) — 



o 



-ifciC 



128 



;i927r 4 ^ 1 + p)C 2 + 0(C 1 )]e 



-27riC 



for 



for c = 0. 



(2.26) 



If we take the limit Re£ — > —00, then the unstable solution V'o (C) decays to zero as a 
power law, according to the expansion (J2.13J) . Thus, the stable solution ^0 (C) nas the 
oscillatory tail given by the representation (J2.21|) with the amplitude factor 

for c 7^ 

(2.27) 



7ikiKi(c) 



for c = 0. 



Similarly, if we take the limit Re^ 



+00, then the stable solution ^0 (0 decays to 



zero, while the unstable solution ^ (C) nas the representation (J2.21j) with the amplitude 
factor given by the negative of expression (|2.27jl . By comparing the other terms on the 
right-hand side of ()2.26|) to the corresponding terms in (|2.24j) . a(c) and p can be uniquely 
determined. 

We now match the leading-order singular behaviour of Vq(p) near p = rtifci, given 
by ()2.25|) . to the formal power series (|2.18|) . Expanding the expressions in ()2.25|) as 
power series gives us 



V (p) 



/fi(c)-a(c) ln(fc?) + 

n=l 

l) n (n + 2)(n + 



0- c 



n 



V 



2n 



E 

n=0 

iiifci. 



(2vr; 



2n 



[(n + 3)5i + p + 



P 



2/i 



for c 7^ 



for c = 0, 



(2.28) 



as p — > ±ik\. These series converge for all \p\ < k\\ in particular, they are valid 
for p — > ±iki, provided \p\ < k\. Hence we can replace (I2.25|) with (|2.28|) in this 
neighbourhood. In (|2.28jh the ellipses stand for coefficients of the expansion of terms 
with a slower growth as p — > ±ik\ which were dropped in (|2.25jl . The discarded terms 
would modify the coefficients of the power series (|2.28|) : however, there are terms which 
would not be affected by these modifications, namely terms with large n. For example, 
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the coefficients proportional to o~(c) and p in (|2.28|) are a factor of n smaller than those 
proportional to Ki(c) and Si; the discarded coefficients would be even smaller. Therefore 
the leading singular behaviour of Vo(p) as p — > ±i&i is determined just by the large- n 
coefficients of the power series (12.28(1 . and hence only the large- n coefficients should be 
matched to the coefficients of the expansion (j2.18(l . This gives the Stokes constant as a 
limit of the coefficients Vm of the series (12.18(1 : 



This formula is used in the next section for numerical computations of the leading-order 
Stokes constant K\{c) for c ^ 0. 

Note that, since (|2.18jl matches (|2.28jl in the limit n — > oo, the formal power series 



diverges for all ( unless Vo(p) converges everywhere (which requires that all the Stokes 
constants be zero). 

Next, we note that as c — * 0, the Stokes constant K\(c) does not tend to Si, its 
value at c = 0. This discontinuity is due to the fact that, as c — > 0, pairs of simple roots 
in the resonant set 1Z C coalesce. (E.g. iki coalesces with ik 2 at 27ri, and so on.) As a 
result, all roots are double and the representation of <pi(C) is discontinuous at c = 0, with 
the power degree r of the prefactor in (j2.23() jumping from r = for c ^ to r = 3 for 
c = 0. In particular, in exceptional models, i.e., discrete models with continuous families 
of stationary kinks (like p. 8(1 . (J1.9j) and (jl.lljl ) the constant Si is a priori zero while the 
limit of Ki(c) as c — > may be nonvanishing. In fact, numerical computations of the 
top limit in ((2.29(1 indicate that the Stokes constant blows up as c — > 0. Renormalisation 
of Ki(c) for small c is, however, a nontrivial asymptotic problem which is beyond the 
scope of our current investigation. 

For c G (ci, 1), where c\ ~ 0.22, the resonant set TZ C contains only one root iki 
and, therefore, there is just one Stokes constant Ki(c), which completely determines 
the convergence of the formal power series for ^o(C)- If ^i(co) = at some point 
c G (ci,l), the stable and unstable solutions i/Jq\C) and V'o^C) coincide to leading 
order. Arguments based on the implicit function theorem (see |TP05| ) reveal a 
heteroclinic bifurcation which occurs on crossing a smooth curve c = c*(h) on the 
(c, /i)-plane, with c*(0) = cq. 

On the other hand, for c G (0, Ci) the resonant set TZ C contains more than one root. 
If Ki(c ) = for some c G (0, Ci), this alone is not sufficient for the convergence of the 
formal power series ^o(C)- The higher-order Stokes constants K 2 (c), K 3 (c), . . . , must 
be introduced and computed from the asymptotic behavior of the power series V (p). 

As we shall show in the next section, the function K\{c) does have zeros in the case 
of the discretizations (Jl. 8(1 and (|1.11(1 . All these zeroes are "safe"; that is, all Cq values 
lie in the interval (ci, 1), so that the higher-order Stokes constants do not have to be 
computed. 




for c 7^ 



for c = 0. 



(2.29) 



Vo(p) also has radius of convergence k\. However, the formal inverse-power series V'o(C) 
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3. Numerical computations of the Stokes constant 

In this section, we report on the numerical computation of the Stokes constants K\ (c) 
for the four different discretizations of the <j) A model (jl.3|) under consideration. Our 
numerical method utilises the expression ()2.29)1 of the Stokes constant in terms of the 
coefficients of the formal power series solution 1)2.18)1 . First, we obtain the recurrence 
relation for the coefficients in the set {f n }^Lo by substituting the power series expansion 
(12.18)1 into the limiting integral equation 1)2.17)1 . and using the convolution formula 

P n *p m = ( " !m | ^ P n+m+l - (3.1) 
[n + m + 1J! 

After that, we compute the asymptotic behavior of these coefficients as n — > oo and 
evaluate the limit (j2.29j) numerically for a fixed value of c ^ 0. 

In order to calculate the Stokes constant for the four models in a uniform way, we 
write a general symmetric homogeneous cubic polynomial Q(u n -i, u n , u n+ i) as 
ill 

Q= ^ ^ ^ y ^ y Q Q,/3,7 u n+a u n+/3 u n+yi (3-2) 
a=— 1 /3=a 7=13 

where a a ^^ are numerical coefficients, with a, (3, 7 G { — 1,0, 1} and a < j3 < 7. The 
symmetry implies that a Qj/ g i7 = a- 7j -p- a and therefore it is sufficient to specify just six 
coefficients. The values of these coefficients for the four nonlinearities in question are 
given in Tabled 



Model 


Ol,l,l 


a 0,0,0 


a 0,l,l 


a 0,0,l 


cJ-1,1,1 


0-1,0,1 


One-site (fTBJ) 





1/2 














Speight- Ward (fTsjl 


1/36 


1/9 


1/12 


1/12 








Bender- Tovbis ()1.9j) 











1/4 








Kevrekidis (fTTT]) 


1/8 











1/8 






Table 1. The coefficients a a ./3 n ~ a f .-f3.- a of the cubic polynomial ()3.2)l for the 

four models under consideration. 

By applying the Borel-Laplace transform ()2.15)) to equation ()2.14)) with Q as in 
([3.2)1 . we obtain the corresponding cubic convolution function Q[V(p)] on the right hand 
side of the integral equation ()2.17j) : 

111 

q [v(p)} = E E E eapy ^ * e ^ PV ^ * e " PV ^- ( 3 - 3 ) 

a=— 1 f3=a 7=/3 

To derive the recurrence formula for the coefficents V2 n in ()2.18j) . it will be more 
convenient to consider the power series expansion which consists of both even and odd 
powers of p: 

v (p) = J2 v ™p n - ( 3 - 4 ) 

n=0 
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We now substitute the series (|3.4|) into (|2.17|) with Q \V~{p)] given by (|3.3|) and use the 
convolution formula ()3.1|) . Equating the coefficients of p n+2 where n — 0, 1, 2, . . ., in the 
resulting equation, we find that 

[n/2] Q 111 ( n /n-i k 

El 2 \ ^ \ ^ \ ^ °a,/3,7 I \ ^ / \ ^ 01 

(2^2)!^ - ct '«=LLL (n + 2)(n+l) \ M ^ 

i=0 V ' a=-l/3=a 7 =/3 V > K ' K i=0 \k=0 



Z. ifi-i Z. ^ 



j=0 \ 1=0 / \m=0 



i\(n — i)\ 



nl 



(3.5) 



where [n/2] is the integer part of n/2 and 0° = 1. Equation (j3.5j) is a recurrence 
relation between the coefficients {w n }^o- Solving equation (|3.5j) with n = 0, we 
get t; = 2-\/i — c 2 . Note that this result is independent of the choice of a Q ,/3, 7 , i.e. 
independent of the model. Letting v\ = and making use of the symmetry of Q, one 
can show by induction that the coefficients of all odd powers in (J3.4j) are zero (as we 
concluded previously on the basis that the outer expansion is odd). 

To prevent overflow or underflow when evaluating the recurrence relation 
numerically, we shall work with the normalised coefficients 

w n = (-l) n kl n v 2n 

so that the Stokes constant ()2.29j) for c ^ is given by 

Ki(c) = lim w n . (3.6) 

n— >oo 

Reformulating ()3.5|) in terms of w n , we use the relation ()3.6|) to compute w n numerically. 
We truncate the sums involving 1/ (2i + 2)!, l/l\ and 1/m! when these factors become 
smaller than 10~ 50 , and evaluate the sums involving the combinatorial factors in two 
halves. In the first, the summation index increases from zero to the halfway point, 
and in the second it decreases from its maximum. This ensures that the combinatorial 
factors are always decreasing from one step to the next so that they can be accurately 
determined recursively. We also truncate these sums when the combinatorial factors fall 
below irr 50 . 

These expedients result in a numerical routine fast enough to allow for evaluation 
of the recurrence relation up to very large n; this is essential given the slow convergence 
of w n to a constant. Matching ()2.18|) to (|2.28j) yields 

v 2n -> (-ly-k^lK^c) + cr(c)/n} as n -» oo; 

therefore, the rate at which w n converges to Ki(c) is of order 1/n: 

Wn = Kl (c) + ^ + *-4 + o(±). (3.7) 
n n z \ n / 

Although the convergence of w n to Ki(c) is extremely slow, we can accelerate the process 
by using ([3.7)1 . Defining 

w n = w n + n(w n - w n -i), 
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2 xlO d 




5 xlO 4 



Figure 2. Convergence of the sequence w n and its accelerated counterpart w n to the 
Stokes constant Ki{c). The dashed line depicts the values of w n and the solid line 
marks the accelerated sequence w n . Left: the Kevrekidis discretization 1)1.11(1 with 
c = 0.5. Right: Speight- Ward discretization 1)1.8(1 with c = 0.005. 



we get 



, \ o-(c) + a(c) ^ ( 1 

W n = KAc) - 1 ' - 1 > +0[ — 



n 2 V n 3 



The convergence of the sequence w n is much faster than that of w n ; see Figure El The 
relative error 

_ ^(c) + cr(c) 1 

can be written as 

E[n) = 

plus terms of order 1/n 4 . This gives an empirical criterion for the termination of the 
process. We continued our computations until E(n) reached a value smaller than 10~ 3 , 
i.e. until the percentage error dropped below 0.1%. For c > 0.5, the value of n to 
which we have to compute in order to achieve this accuracy is less than 100, increasing 
for smaller values of c to approximately 5 000 for c = 0.005. Consequently, the above 
numerical algorithm is not suited to the study of the c — >• limit, and would have to be 
modified for that purpose. 

Figure El displays the Stokes constant computed using the above numerical 
procedure, for the four models of Table HJ We see that the Stokes constant Ki(c) 
vanishes almost nowhere in c ^ in all four models. There are, however, several 
isolated zeros: Ki(c) = for Co ~ 0.45 in the case of the Speight- Ward nonlinearity 
(|1.8|) and for cq ~ 0.37, 0.63 and 0.83 in the case of the Kevrekidis discretization (jl.lljl . 
Importantly, all of these lie in the region (ci, 1) where the resonance set (|2.2()|) consists 
of only one value, p\ = iki. (Here C\ ~ 0.22.) Therefore, there is a sliding kink in the 
h — >• limit for each of these isolated values of velocity. Furthermore, strong parallels 
between our current setting and that of solitons of the fifth-order KdV equation |TP05j 
suggest that sliding kinks should exist along a curve on the (c, h) plane emanating from 
each of the points (cq, 0). In other words, we conjecture that there is a radiationless kink 
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Figure 3. The Stokes constant K\ as a function of the kink's speed c for the four 
discretizations of the </> 4 model. Clockwise from the top left: one-site ljl.6|l : Bender- 
Tovbis (TJj!; Kevrekidis lfl~TT|) : Speight- Ward JTHJl 



travelling with a certain particular speed c*(/i) for each h in the case of the Speight- 
Ward nonlinearity, and that there are three such velocities (for each h) in the case of 
the Kevrekidis model. For small h, c*(h) should be close to the above values Cq. 

In order to verify the existence of kinks sliding at these isolated velocities by an 
independent method, we solved the differential advance-delay equation ()2.2|) numerically. 
The infinite line was approximated by an interval of length 2L = 200, with the 
antiperiodic boundary conditions <f)(—L) = — 4>{L). We utilised Newton's iteration with 
an eighth-order finite-difference approximation of the second derivative; the step size 
was chosen to be h/10. The continuum solution (jl.2j) was used as an initial guess. 

If we find a solution to the advance-delay equation with (p(z) decaying to a 
constant for large positive and negative z, then we regard this solution as (a numerical 
approximation to) a radiationless travelling kink. We were able to tune c for a fixed 
value of h so that the radiation was reduced to the order of 10~ 12 , whereupon the finite 
accuracy of our numerical scheme prevented any further reduction. To make sure that 
the radiation does vanish rather than reaching a local minimum but remaining nonzero, 
we plot the average magnitude of the radiation near the ends of the interval as a function 
of c, for fixed h. This is defined as the average of \_4>{z) — 0] over the last 20 units of 
the interval, where is the average value of <fr(z) over these last 20 units. The results are 
shown in figure EJ Note the straight-line behavior of the graphs near the isolated values 
of c; this indicates that the coefficient of the sinusoid superimposed over the kink's flat 
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Figure 4. The numerical evidence for the disappearance of radiation at isolated 
values of c = c* in Speight and Ward's model (top left panel) and Kevrekidis' model 
(three other panels). The magnitude of the oscillatory tails, as defined in the text, is 
plotted as a function of c for a fixed value of h. The minimum radiation is attained at 
the value c*(/i) which is found numerically. (This c*(/i) is of course slightly different 
from the value c*(0) for which the Stokes constant vanishes.) 



asymptote crosses through zero (rather than attaining a small but nonzero minimum). 
The supression of radiation at the isolated points is thereby verified. 

Finally, the last question that we would like to address here is whether the intensity 
of the radiation from the moving discrete kink depends on the type of discretization. 
More specifically, we would like to know whether the choice of one of the exceptional 
discretizations (which, by definition, support translationally invariant stationary kinks) 
serves to reduce the radiation from the moving kinks. Speight and Ward have already 
given an affirmative answer for their exceptional discretization; here we consider the 
one-parameter nonlinearity 

Q = ^y^n + ~ u 2 n {u n+1 + Mn _0, (3.8) 

which interpolates between the one-site nonlinearity (jl.6|) (for which \i = 0) and the 
exceptional discretization (jl.9|) of Bender and Tovbis (for which \l — 1). Figure El shows 
the Stokes constant for the model (jll8|1 . as /i changes from to 1 for fixed values of 
c. The Stokes constant is indeed seen to be drastically reduced as fi approaches 1 - 
that is, in the limit of the exceptional discretization. (It nonetheless remains nonzero, 
of course, unless c = 0.) 
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Figure 5. The Stokes constant as a function of fi in the model (|3.8|) . Note the 
logarithmic scale of the vertical axis. 

4. Concluding remarks and conclusions 

In this paper we have investigated the existence of sliding kinks — i.e. discrete kinks 
travelling at a constant velocity over a flat background, without emitting any radiation 
- in four discrete versions of the quartic-coupling theory. One of these models is the 
most common, one-site, discretization. As the overwhelming majority of discrete <fr A - 
equations, it supports travelling kinks, but these kinks do radiate and decelerate as a 
result. The other three discretizations we considered are all exceptional in the sense 
that they all support one-parameter continuous families of stationary kinks where the 
free parameter defines the position of the kink relative to the lattice. This property is 
clearly nongeneric; the translation invariance of the continuous 4 -theory is broken by 
the discretization and hence in generic disretizations kinks may only be centred at a 
site or midway between two sites. Since the nonexistence of "translationally-invariant" 
and of sliding kinks in the generic models can be ascribed to similar factors, viz., the 
breaking of the translation and Lorentz invariances, it was hoped that the exceptional 
discretizations might turn out to be equally exceptional from the point of view of sliding 
kinks. Our approach was based on the computation of the Stokes constants associated 
with the putative sliding kink in a given equation. 

The main conclusion of our work is that the sliding kinks do exist in the discrete 
<f) A theories, but only with special, isolated, velocities (which of course depend on h). 
There is one such velocity in the exceptional model of Speight and Ward, and three 
different sliding velocities in the discretization of Kevrekidis. It is natural to expect 
that the sliding kinks should play the role of attractors similarly to the fronts moving 
with "stable velocities" in dissipative systems; that is, radiating travelling kinks should 
evolve into kinks travelling with the sliding velocities — if there are such velocities in 
the system. Not every discretization supports sliding kinks, of course; in particular, no 
sliding velocities arise for the generic, one-site, nonlinearity and even for the exceptional 
discretization of Bender and Tovbis. 

One natural way of trying to construct the sliding kinks is via power series 
expansions in powers of c 2 ; for the exceptional discretizations, this construction can 
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be carried out to any order. This approach was pursued in the recent work of Ablowitz 
and Musslimani |AM03j . Our results indicate, however, that these power series will 
not converge and exponentially small terms (terms lying beyond all orders of the power 
expansion) emerge because of the singular behaviour of the Stokes constant Ki(c) as 
c — > 0. Detailed studies of this singular limit will be presented elsewhere. 

The exceptional discretizations have richer underlying symmetries than generic 
nonlinearities but the "translation invariance" of the stationary kink alone does not 
automatically guarantee the existence of the sliding velocities. The exact relation 
between the "translational invariance" and mobility of kinks is still to be clarified; at 
this stage it is worth mentioning that the Stokes constant associated with (and hence 
the intensity of radiation from) a moving kink is several orders of magnitude smaller in 
exceptional models than in generic discretizations. 

Finally, it is instructive to point out some parallels with an earlier work of Flach, 
Zolotaryuk and Kladko |FZK99j who also studied the phenomenon of kink sliding in 
Klein-Gordon lattices. In the scheme of |FZK99j , one postulates an analytic expression 
for the sliding kink, u n {t) = <p{n — ct — s), with some explicit function 4>(z), and then 
reconstructs the Klein-Gordon nonlinearity for which this is an exact solution. Our 
present conclusions are in agreement with the results of these authors who observed 
that for a given h, the kink may only slide at a particular, isolated, velocity. The two 
approaches, ours and that of [FZK99 , are reciprocal; while we examine the existence 
of sliding kinks for particular discretizations of the 4 -theory, with fixed parameters 
independent of the kink's velocity, in the "inverse method" of |FZK99j one assumes 
an explicit solution of a particular form but does not have any control over the 
resulting nonlinearities. Consequently, the discrete Klein-Gordon models generated by 
the "inverse method" are not discretizations of the 4 -theory and do include explicit 
dependence on the velocity of the sliding kink. 
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